Gravitating Monopole-Antimonopole Systems at 
Large Scalar Coupling 

Jutta Kunz, Ulrike Neemann, Yasha Shnir 

Institut fiir Physik, Universitat Oldenburg, Postfach 2503 
D-26111 Oldenburg, Germany 

February 2, 2008 

Abstract 

We discuss static axially symmetric solutions of SU(2) Einstein- Yang- Mills- 
Higgs theory for large scalar coupling A. These regular asymptotically flat solu- 
tions represent monopole-antimonopole chain and vortex ring solutions, as well as 
new configurations, present only for larger values of A. When gravity is coupled 
to the Yang-Mills-Higgs system, branches of gravitating solutions emerge from the 
flat-space solutions, and extend up to critical values of the gravitational coupling 
constant. For small scalar coupling only two branches of gravitating solutions 
exist, where the second branch connects to a generalized Bartnik-McKinnon solu- 
tion. For large scalar coupling, however, a plethora of gravitating branches can be 
present and indicate the emergence of new flat-space branches. 

PACS numbers: 14.80.Hv,11.15Kc 
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1 Introduction 



The 't Hooft-Polyakov monopole pQ represents but the simplest of a rich variety of regu- 
lar classical non-perturbative finite energy solutions of SU(2) Yang-Mills-Higgs (YMH) 
theory. Besides this spherically symmetric solution, axially symmetric multimonopoles 
[21 El HI E], monopole-antimonopole pairs (MAPs) [HI [7J E], monopole-antimonopole 
chains (MACs) [9] , and vortex ring solutions [HI [10] are known, as well as multimonopole 
solutions with only discrete symmetries [TT] . 

In the Bogomornyi-Prasad-Sommerfield (BPS) limit of vanishing Higgs potential, 
the repulsive and attractive forces between two monopoles exactly compensate for any 
separation, thus BPS monopoles experience no net interaction [12]. As the scalar field 
becomes massive, the fine balance of forces between the monopoles is broken, since the 
attractive Yukawa interaction becomes short-ranged. Consequently non-BPS monopoles 
experience repulsion [5j. This repulsion can be overcome by the inclusion of gravity, 
which allows for bound multimonopoles [T3] . 

In monopole and multimonopole solutions, the nodes of the Higgs field are associated 
with the location of the magnetic charges. This also holds for monopole-antimonopole 
pair and chain solutions, where m monopoles and antimonopoles are located symmetri- 
cally on the symmetry axis in alternating order [6J [7J El [9]. When the magnetic charge 
n of the individual monopoles and antimonopoles becomes large, the balance between 
the repulsive and attractive interactions of the monopoles and antimonopoles is shifted. 
For vanishing Higgs potential, the structure of the solutions then changes completely, 
and the zeros of the Higgs field form one or more rings, centered around the symmetry 
axis. Thus vortex ring solutions arise |9J. 

When the scalar coupling A is small, the transition from MAPs and MACs to vortex 
ring solutions appears for magnetic charge n = 3 [9]. For larger scalar coupling A, 
however, bifurcations arise, where new pairs of solutions appear [10]. Thus MAP or 
MAC solutions and vortex ring solutions may coexist beyond some critical value of 
A. Moreover new configurations with mixed node structure arise, possessing isolated 
nodes on the symmetry axis as well as vortex rings [10]. The strength of the Higgs self- 
interaction is clearly critical for the existence and the properties of these new equilibrium 
configurations. 

We here address such monopole-antimonopole systems in flat-space, reviewing and 
supplementing previous results [10]. In particular, we obtain new sets of solutions as- 
sociated with further bifurcations at larger values of the scalar coupling, both in the 
topologically trivial and the non-trivial sector. While the mass of these configurations 
increases with the scalar coupling, it is largely determined by their internal structure. We 
confirm our previous conjecture, that for large scalar coupling, MAC solutions appear 
to be the energetically most favourable configurations. 

The balance of forces changes again, when gravity is coupled to YMH theory. This 
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makes the investigation of self-gravitating MAPs, MACs and vortex ring solutions of 
Einstein- Yang-Mills-Higgs (EYMH) theory very interesting. Depending on the strength 
of gravity, the structure of the solutions may change again [Hj. Also new solutions may 
arise, which do not possess a flat-space limit [T5| [T6| [T7] . But gravity cannot become 
too strong, for such regular equilibrium configurations to exist, since beyond a certain 
gravitational coupling strength, a horizon is expected to form |15j . 

For gravitating monopoles and multimonopoles a degenerate horizon forms indeed 
at a critical value of the gravitational coupling [151 US 13] . For monopole-antimonopole 
systems, however, no horizon is seen to form. Instead the flat-space branch of regular 
graviting monopoles bifurcates with a second branch at a maximal value of the grav- 
itational coupling. This second branch then extends backwards and, in the limit of 
vanishing coupling constant a, it becomes linked to a generalized Bartnik-McKinnon 
(BM) solution of Einstein- Yang-Mills (EYM) theory [T9| 1201 [21] . 

Anticipating new phenomena [181 121] , we here address the effect of gravity on mono- 
pole-antimonopole systems at large values of the scalar coupling. Indeed, as A becomes 
large, we observe a rich pattern of bifurcating branches of solutions arising at finite 
gravitational coupling strength, which encompasses all possible combinations of isolated 
nodes and vortex rings. Moreover, certain bifurcations of the gravitating branches of 
solutions turn out to be precursors of the bifurcations of the flat-space solutions, and, in 
the limiting case of infinite A, particular bifurcations become linked to the corresponding 
generalized BM solutions of EYM theory. 

In section II we present the action of EYMH theory, the axially symmetric ansatz 
and the boundary conditions. We then discuss the properties of gravitating monopole- 
antimonopole paris, chains and vortex rings at large A in section III, focussing on 
monopole-antimonopole systems with m = 2, 3, 4 and n = 3. We give our conclusions in 
section IV. 



2 SU(2) EYMH Action and Axially Symmetric Ansatz 
2.1 Action 

We consider the SU(2) Einstein- Yang-Mills-Higgs theory with action 
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with curvature scalar R, su(2) field strength tensor 



d^Ay — dyA^ + ie[A^, A v \ , 



(2) 
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gauge potential = A'j 1 T a /2, and covariant derivative of the Higgs field $ = $ a r a in 
the adjoint representation 

D^ = d^ + ie[A^]. (3) 

Here G and e denote the gravitational and gauge coupling constants, respectively, 77 
denotes the vacuum expectation value of the Higgs field, and A represents the strength 
of the scalar coupling. 

Under SU{2) gauge transformations U, the gauge potentials transform as 

A^ = UA^ + -(d,U)U^ , (4) 

e 

and the Higgs field transforms as 

$' = U$U ] . (5) 

The nonzero vacuum expectation value of the Higgs field breaks the non-Abelian 
SU(2) gauge symmetry to the Abelian U(l) symmetry. The particle spectrum of the 
theory then consists of a massless photon, two massive vector bosons of mass M v = erj, 
and a massive scalar field M s = y/2Ar}. In the BPS limit the scalar field also becomes 
massless, since A = 0, i.e., the Higgs potential vanishes. 

Variation of the action ([!]) with respect to the metric g^ u leads to the Einstein 
equations 

Gp,v = Rnu ~ -jdixvR = SttGT^ (6) 

with stress-energy tensor 

rp TO 



QgflV 

= 2 Tr {F m F v ^ - l -g^F aP F aP ) 

+ Tr {\d^D u § - ^g^D^D"*) - ^^Tr($ 2 - r, 2 ) 2 . (7) 

Variation with respect to the gauge field A^ and the Higgs field $ leads to the matter 
field equations, 

1 ^(v^n-k^.^Ho, (8) 



1 



respectively. 



D M ( x /^ J D^) + A($ 2 -r ? i )$ = , (9) 
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2.2 Static axially symmetric Ansatz 

To obtain gravitating static axially symmetric solutions, we employ isotropic coordinates 
[20, [131 In terms of the spherical coordinates r, 9 and ip the isotropic metric reads 

, o „ _ n m , o mr 2 Zr 2 sin 2 # , 9 , 
rf s 2 = _/dt 2 + _rfr 2 + — d0 2 + V , (10) 

where the metric functions /, m and I are functions of the coordinates r and 0, only. 
The z-axis (0 = 0, ir) represents the symmetry axis. Regularity on this axis requires 

m|e = o l7 r = l\e=o,n ■ (11) 

We take a purely magnetic gauge field, A = 0, and parametrize the gauge potential 
and the Higgs field by the Ansatz [H] 



A^dxP = \—dr + (1 - K 2 )de^ Tlp 



in) 



\ r J 2e 

((n,m) (n,m) \ 

K 3 ^— + K 4 ^—\d^, (12) 

$ = v (W r ^ m ) + $ 2 r, (n ' m) ) . (13) 

where the su(2) matrices Tr U ' m \ Tg n ' m \ and l~!p ^ £1X6 defined as products of the spatial 
unit vectors 

_ ( s i n ( m #) cos(nif), sm(m6) sin(n^), cos(m8)) , 
eg" 1 ' = (cos(m6) cos(n^), cos(m9) sin(n^), — sin(m6>)) , 
e^ n) = (-sin(n^),cos(n^),0) , (14) 

with the Pauli matrices r a = (r x , r y , t z ), i.e. 

T (n,m) = sin ( m 0) T M + C os(m0)r 2 , 
T (n,m) _ coa ( m Q'j T j.ri) _ S in(m^)r x , 
r£ n) = - sm(rap)r x + cos(rnp)r y , 

with Tp = cos(nip)r x + sin(ray9)r y . For m — 2, n = 1 the Ansatz corresponds to the 
one for the monopole-antimonopole pair solutions [3, 0, US], while for m = 1, n > 1 it 
corresponds to the Ansatz for axially symmetric multimonopoles [3j El [13]. The four 
gauge field functions Ki and two Higgs field functions $j depend on the coordinates r 
and only. All profile functions are even or odd w.r.t. reflection symmetry, z — > —2. 
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The gauge transformation 

U = exp{zT(r,9)r^/2} (15) 

leaves the Ansatz form-invariant [22]. To construct regular solutions we have to fix the 
gauge [5]. Here we impose the gauge condition [9] 

rd T K x - d e K 2 = . (16) 

With this Ansatz the equations of motion reduce to a set of 9 coupled partial dif- 
ferential equations, to be solved numerically subject to the set of boundary conditions, 
discussed below. 

2.3 Boundary conditions 

To obtain globally regular asymptotically flat solutions with the proper symmetries, we 
must impose appropriate boundary conditions [T51 19"]. 
Boundary conditions at the origin 

Regularity of the solutions at the origin (r = 0) requires for the metric functions the 
boundary conditions 

d r f(r,6)\ r=0 = d r m(r,9)\ r=0 = d r l(r,9)\ r=0 = , (17) 

whereas the gauge field functions Ki satisfy 

K x (0, 0) = K 3 (0, 9) = K 4 (0, 9) = , K 2 (0, 9) = 1 , (18) 

and the Higgs field functions satisfy 

sin(m0)$i(O, 9) + cos(m#)$ 2 (0, 9) = , (19) 

d r [cos(m0)$i(r,0) - sin(m0)$ 2 (r, 0)] \ r=Q = , (20) 

i.e. $,(0,0) =0, <9 r $,(O,0) = 0. 
Boundary conditions at infinity 

Asymptotic flatness imposes on the metric functions of the solutions at infinity (r = 
oo) the boundary conditions 

/ — I, m — 1, Z— >1 . (21) 

Considering the gauge field at infinity, we require that solutions in the vacuum sector 
Q = 0, where m = 2k, tend to a gauge transformed trivial solution, 

$ ^ v Ur z tf , — > Up^lfl , 
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and that solutions in the sector with topological charge n, where m = 2k + 1, tend to 
$ — Eft , A M — > £Mg£W + ^(^C/)C/t , 



T W T (i.») 

' to , ^ . ^ ' a 



= , = -^-dfl - n sin 6>-^— dy? 



where 

4>(1,«) _ r)r (V0 4(i,n) r ; 7 ./* _ 

2e 2e 

is the asymptotic solution of a charge n multimonopole, and U = exp{—ik6r^}, both 
for even and odd m. 

In terms of the functions K\ — K4, $1, $2 these boundary conditions read 

K x — > , K 2 — > 1 - m , (22) 

cos 6> - cos (ra0) l-cos(m0) 

a 3 — ► m odd , A 3 — > m even , (23) 

sin sin 9 

sin 

$! — > 1 , $ 2 ^0. (25) 
Boundary conditions along the symmetry axis 

The boundary conditions along the z-axis (6 = and 9 = 71") are determined by the 
symmetries. The metric functions satisfy along the axis 

d e f = d e m = d e l = , (26) 

whereas the matter field functions satisfy 

Kt = K 3 = $ 2 = , 9^2 = d e K 4 = <%$i = . (27) 

2.4 Mass and charge 

Let us introduce the dimensionless coordinate x and the dimensionless coupling constant 
a P], 

ea 

x = =r , a = V47rGr7 . (28) 
yAirG 

The limit a — ► can be approached in two different ways: G — » while the Higgs 
vacuum expectation value 77 remains finite (flat-space limit), and 77 — > while Newton's 
constant G remains finite. Corresponding to these two limits two branches of solutions 
may exist. 

The dimensionless mass M of the solutions is obtained from the asymptotic expansion 
of the metric function /, 

M = -L lim x 2 d x f . (29) 
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The asymptotic expansion of the gauge fields yields the dimensionless magnetic charge 
P of the solutions, 

P = | [1 " (-1)"] , (30) 

i.e., solutions in the topologically trivial sector have no charge, P = 0, whereas solutions 
in the non-trivial sectors have charge P = n. 

3 Numerical Results 

We have constructed numerically regular gravitating monopole-antimonopole systems, 
including monopole-antimonopole chain and vortex ring solutions. We here focus on 
solutions with m = 2, 3 and 4 and n = 3. In particular, we illustrate the plethora of 
solutions appearing, when the scalar coupling A is sufficiently large. At the same time we 
study the dependence of the properties of these solutions on the gravitational coupling 
strength a. 

We first briefly address the numerical procedure. Then we present our results for 
monopole-antimonole chains, vortex rings and configurations with more complicated 
structure, beginning with the m = 2, n = 3 systems. 

3.1 Numerical procedure 

To construct solutions subject to the above boundary conditions, we map the semi- 
infinite range of the radial coordinate r onto the closed unit interval of the new com- 
pactified radial variable x G [0:1], 

T 

X = , 

1 + r 

i.e., the partial derivative with respect to the radial coordinate changes according to 

d r -> (1 - Xfdx . 

The numerical calculations are then performed with help of the package FIDISOL, based 
on the Newton- Raphson iterative procedure [23] . Technically, the system of non-linear el- 
liptic partial differential equations of second order is first discretized on a non-equidistant 
grid in x and 9. Typical grids used have sizes 60 x 40 covering the integration domain 
< x < 1 and < 6 < tt/2. The numerical errors are typically of order of 10~ 4 . 

3.2 MAPs and vortex rings: rn = 2. n = 3 
3.2.1 Flat space solutions 

Monopole-antimonopole pairs (MAPs) reside in the topologically trivial sector. These 
solutions represent unstable equilibrium configurations and can be thought of as con- 
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sisting of two solitons of individual charges ±n. The position of each constituent can 
be identified according to the location of the nodes of the Higgs field. The separation 
between the poles depends on the subtle interplay of forces and thus on the strength of 
the scalar coupling A [8]. 

The MAPs with n = 1 exist for arbitrary values of A. They have two isolated nodes 
on the symmetry axis, indicating the position of the poles, and their energy density 
possesses two local maxima, associated with the nodes of the Higgs field [TB] . 

As the charge of the poles increases, the interaction between the nonabelian matter 
fields becomes stronger than in the case of unit charge constituents, leading to new types 
of solutions. Indeed, for small and vanishing values of A, MAPs with n > 2 no longer 
exist. Instead solutions of vortex ring type emerge, where the Higgs field vanishes on a 
ring centred around the symmetry axis in the xy plane [H] . These flat-space vortex ring 
solutions form the fundamental A-branches, which persist for arbitrary values of A. As 
the scalar coupling strength increases, however, bifurcations are observed, where pairs 
of additional flat-space solutions emerge at critical values of A |10j . 

For n = 3, in particular, a pair of MAP solutions arises at the critical value A* = 
1.382. The A-dependence of these flat-space solutions is demonstrated in Fig. [Tj, where 
we show their mass and the location of the nodes of their Higgs field configurations. 
At the bifurcation, the mass of the new MAP solutions is higher than the mass of the 
vortex ring solution. For large values of A, however, the MAP solution on the lower mass 
branch becomes energetically favoured. Note, that the node structure of the solutions 
on the upper mass branch changes from MAP to vortex ring, as A increases, but the 
radius of the ring remains small. 
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(b) ° 2 4 x'« 6 

Figure 1 : (a) The mass M of the flat-space solutions on the fundamental A-branch as well as 
on the two additional A-branches for m = 2, n = 3 solutions, (b) The location of the isolated 
nodes z and the radius p of the vortex ring of the Higgs field for the same set of solutions. (See 
Table 1 for the notation of the node and vortex ring configurations of the Higgs field.) 
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MA 






A monopole-antimonopole pair with a 
monopole above and an antimonopole below 
the xy-plane on the symmetry axis. 






V 






A vortex ring in the xy-plane. 



Table 1: Configurations of nodes and vortex rings of the Higgs field for m = 2, n = 3 solutions. 
3.2.2 Gravitating solutions 

When gravity is coupled, branches of gravitating solutions arise from these flat-space 
configurations. For sufficiently large gravitational coupling a, however, a gravitational 
instability should develop. Therefore regular gravitating solutions will exist only up to 
a maximal value of a. 

The critical behaviour of these a-branches, emerging from flat-space solutions, is 
elucidated by numerical investigation. As such an a-branch approaches its maximal 
value of a, it merges with a second a-branch of regular gravitating solutions, where 
ttmax is the endpoint of both branches. The nature of the second branch, however, 
depends on the nature of the first branch and on the value of A. 

Fundamental solutions 

When considering sets of solutions, emerging from the fundamental flat-space A- 
branch, we always find two a-branches, where the solutions on the second a-branches 
shrink to zero size in the limit a — > 0, while their mass diverges. By scaling the coordi- 
nates and the Higgs field of the solutions via 

x = - , $ = ct$ , (31) 

a 

one realizes that the second a-branches connect to limiting EYM solutions, when a — > 0, 
solutions with finite scaled size and finite scaled mass M jT6] , 

M = olM . (32) 

Indeed, after the scaling, the field equations do not depend on a. Instead a appears 
in the asymptotic boundary conditions of the Higgs field, |$| — > a. The Higgs field then 
becomes trivial on the second a-branches as a — > 0, and the solutions thus approach 
EYM solutions. Depending on n, the limiting solutions correspond to the lowest BM 
or a generalized BM solution, i.e., the solutions on these second a-branches possess no 
flat-space limit. 
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The two a-branches are illustrated in Fig. [2] for solutions with m = 2, n = 3 
coupling strength A = 0, A = 2.0 and A = 4.5. The figure exhibits the mass and the 
scaled mass of these solutions and the location of their nodes. Note, that they always 
correspond to vortex ring solutions. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

(a) a (b) a 

Figure 2: (a) The mass M and the scaled mass aM of m = 2, n = 3 solutions, emerging from 
the fundamental flat-space A-branch versus a. (b) The radius p of the vortex ring of the Higgs 
field for the same set of solutions. 

New solutions: intermediate A 

The second a-branch always connects to an EYM solution, when the first a-branch 
emerges from a fundamental A-branch solution. But for the a-branches of the new flat- 
space solutions, which exist only beyond a critical value of A, one may anticipate a 
different pattern, because these solutions always come in pairs. 

In the presence of gravity, each of the solutions of such a pair forms an a-branch, 
extending up to a maximal value of a. But the critical behaviour at a max now depends 
on the strength of the scalar coupling A. In the simplest case, observed for not too large 
values of A, the two a-branches simply merge and end at a max . Indeed, this critical 
behaviour is different from the above, since both branches merging at a max have flat- 
space limits. 

This is illustrated in Fig. [3J where we exhibit the mass and the location of the nodes 
for m = 2, n = 3 configurations at A = 2.0, A = 2.645, A = 3.645, A = 4.5, A = 5.25 
and A = 12.5 versus the coupling constant a. Evidently, for the smaller A values, a first 
a-branch of solutions originating from the lower mass flat-space solution, extends up to 
a maximal value a max , where it merges with a second a-branch, originating from the 
upper mass flat-space solution. The solutions on these branches correspond to MAPs. 

However, as A increases, a new phenomenon arises on the second a-branch, as it 
evolves back from a max towards the flat-space solution. As seen in Fig. [3], at a critical 
value of A (A ~ 4.5) a bifurcation appears on the second a-branch, giving rise to two 
additional a-branches of solutions. 
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— X=2.0: MA 

— X=2.645: MA 

— - A=3.645: MA 

A=4.5: MA,V 
A.=5.25:MAV 
>.= 12.5:MA.V 



0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 0.3 

(a) a (b) a 

Figure 3: (a) The mass M of m = 2, n = 3 solutions, related to the new flat-space A-branches 
versus a for several values of the scalar coupling A. (b) The location of the isolated nodes z 
and radius p of the vortex ring of the Higgs field for the same set of solutions. (See Table 1 
for the notation of the node and vortex ring configurations of the Higgs field.) 

These 3rd and 4th a-branches exist beyond a minimal value of a, cVm> and below 
a respective maximal value of a, where they merge with the first and second a-branch, 
respectively. The range a m i n < amL, where these additional branches exist, increases 
with increasing A. In Fig. [3]the additional branches manifest in the mass of the A = 12.5 
solutions as a higher mass spike. Note, that in this range of A vortex ring solutions arise 
on a part of the second a-branch. 

The emergence and evolution of these new a-branches is illustrated further in Fig. 0], 
where we exhibit the values of the metric functions / and I at the origin, f(0) and 1(0). 
In particular, Z(0) clearly exhibits the various branches and critical points and therefore 
represents an instructive means to clarify the pattern of solutions. 

To complement the previous discussion of the new solutions, let us consider the A- 
dependence of the solutions for fixed values of a. In Fig. [5] we exhibit the mass, the 
location of the nodes, and the values of the metric functions at the origin versus A for 
m = 2, n = 3 configurations at a = 0.2, a = 0.21, a = 0.22, a = 0.23 and a = 0.24. 
One clearly recognizes the new solutions, present in a certain region of parameter space. 

Since both the scalar interaction and gravity are attractive, one can expect a similar 
effect concerning the changing structure of the nodes, when A is increased at fixed a as 
when a is increased at fixed A. For the new flat-space upper mass A-branch, a transition 
from MAPs to vortex rings appears at A = 3.941 [TO] . 

As the gravitational coupling constant a increases, the pattern of the nodes of the 
solutions changes significantly as well, and there are no solutions with vortex rings 
(except for those connected to the fundamental A-branch), when a > 0.215 for any 
value of A. This pattern is illustrated in Figure El where the location of the nodes, both 
isolated zeros and vortex rings, is shown versus the scalar coupling for several values of 



12 




0.3 - 



0.2 0.1 

0.18 0.2 0.22 0.24 



0.05 0.1 0.15 0.2 0.25 

(a) a 



1.2 
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12.5 



0.6 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 

(b) 



0.7 



Figure 4: The values of metric functions at the origin, /(0) (a) and 1(0) (b) of m = 2, n = 3 
solutions, related to the new flat-space A-branches versus a for several intermediate values of 
the scalar coupling A. 



a. Evidently, the coupling to gravity restricts the range of possible values of A, where 
these solutions can exist, and this restriction is the stronger, the larger a, shrinking the 
domain of existence to zero at a critical value of a. 
New solutions: very large A 

As A increases further, the new a-branches extend further backwards. Along these 
branches the mass of the solutions increases further as a decreases. At the same time 
the locations of the nodes of the Higgs field move towards the origin, as illustrated in 
Fig. [6] for very large values of the scalar coupling A. 

With help of Fig. [6] we can now address the limit A —>■ oo for these new a-branches. 
As A increases without bound, the minimal value a m i n , beyond which the new branches 
exist, decreases towards zero. (Note, that a min ~ 1/vA for very large A.) At the same 
time the scaled mass of the critical solutions at a m i n decreases towards the mass of the 
generalized BM solution with n = 3 [20], while their mass itself diverges. Likewise, the 
values of their metric functions at the origin approach those of the n = 3 generalized 
BM solution [21]. Their nodes finally move continuously towards the origin, and the 
solutions shrink to zero size as A — > oo. Scaling the coordinates and Higgs field again 
reveals the generalized BM solution as limiting solution of finite scaled mass and finite 
scaled size [TBJ HZ] . 

This limiting behavior for the critical solutions on the new a-branches is similar to 
that observed for the solutions on the second a-branches of the fundamental solutions. 
In both cases the respective generalized BM solution is approached. 

One can understand the reason for this behaviour here, by noting that, as the scalar 
field becomes infinitely heavy, it decouples in these solutions and, in the double limit of 
infinite A and zero n, it does not affect the dynamics of the gauge sector anymore. In 
that sense the EYMH system is getting actually truncated here to an EYM system. 
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Figure 5: (a) The mass M of m = 2, n = 3 solutions, related to the new flat-space A-branches 
versus A for several values of the coupling a. (b) The location of the isolated nodes z and the 
radius p of the vortex ring of the Higgs field , and the values of metric functions at the origin, 
/(0) (c) and 1(0) (d) for the same set of solutions. (See Table 1 for the notation of the node 
and vortex ring configurations of the Higgs field.) 



3.3 Topologically nontrivial sector: m = 3, n = 3 
3.3.1 Flat space solutions 

Let us now turn to solutions in the topologically nontrivial sector with m = 3, n = 3. In 
these solutions a triply charged monopole resides at the origin, providing the topological 
charge of the solutions. Again we first address the flat-space solutions, reviewing and 
supplementing previous results [TO] . 

For A — > 0, the m = 3, n = 3 solution possesses a triply charged monopole at the 
origin and two oppositely oriented vortex rings located symmetrically above and below 
the xy-plane For larger values of A one expects a bifurcation to occur, where new 
branches of solutions appear, which possess a different node structure, with three isolated 
nodes on the symmetry axis. For high values of A such a monopole-antimonopole chain 
(MAC) should then represent the energetically most favourable solution. 
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Figure 6: (a) The scaled mass aM of m = 2, n = 3 solutions, related to the new flat-space 
A-branches versus a for very large values of the scalar coupling A. (b) The location of the 
isolated nodes z and the radius p of the vortex ring of the Higgs field, the value of metric 
functions at the origin, 1(0) (c) for the same set of solutions, (d) The minimal value a m in of 
the new a-branches of m = 2, n = 3 solutions versus A. (See Table 1 for the notation of the 
node and vortex ring configurations of the Higgs field.) 



Indeed, at A Cl = 0.8099 two new branches of solutions appear which possess the 
node structure of MACs [9J. But as seen in Fig. the vortex ring solutions on the first 
(fundamental) A-branch turn already into MACs before A Cl is reached: As A increases, 
two nodes emerge from the origin and separate from each other along the z-axis. The 
solutions then possess three nodes on the symmetry axis and two vortex rings, located 
symmetrically above and below the xy-plane. As A increases further, the new nodes 
move further apart, while the vortex rings shrink to zero size, merging with the new 
nodes on the symmetry axis. 

Thus beyond a certain value of A, which is smaller than the first critical value A Cl , no 
vortex ring solutions are present any longer [9j. The fundamental branch and the new 
upper branch then merge at A C2 = 0.8194, while the new lower branch continues to large 
values of A, representing the fundamental branch there. This is seen in Fig. where 
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Figure 7: (a) The mass M of the flat-space solutions on the fundamental A-branch as well 
as on the additional A-branches for m = 3, n = 3 solutions, (b) The location of the isolated 
nodes z and vortex rings (p, z) of the Higgs field for the same set of solutions. (See Table 2 for 
the notation of the node and vortex ring configurations of the Higgs field.) 



the mass and the node structure of the m = 3, n = 3 A-branches are exhibited, with 
emphasis on the three branches of MACs present in the small range of A, A Cl < A < A C2 . 

Interestingly, a further bifurcation arises at A C3 = 2.3436, which was missed before. 
At A C3 two new branches of vortex ring solutions appear, whose single ring is located 
in the xy-plenae. Thus beyond A C3 , vortex ring solutions are present again, but possess 
higher masses than the fundamental MAC solutions. These new A-branches are also 
exhibited in Fig. 

3.3.2 Gravitating solutions 

When gravity is coupled, branches of gravitating solutions arise from these m = 3, 
n = 3 flat-space configurations. Depending on the value of A, two or more a-branches 
are present, whose nature is considered in detail in the following. 

The a-dependence of the A < A Cl solutions is simple. A first a-branch emerges from 
the flat-space solution and merges at some a max with a second a-branch, which connects 
to the generalized BM solution in the limit a — > 0. Indeed, for small values of a on the 
second branch, the solutions may be thought of as composed of a scaled n = 3 generalized 
BM solution in an inner region and a flat-space charge-3 monopole configuration in an 
outer region [T7]. In the limit a — > 0, the mass of the solutions diverges, while their 
scaled mass approaches the mass of the generalized n = 3 BM solution. 

Solutions in the range A Cl < A < A C2 

While for A < A Cl only a single flat-space solution is present, there are three flat- 
space solutions in the range A Cl < A < A C2 . This complicates the picture considerably. 
We exhibit in Fig. [H] the node structure of the a-branches, associated with all three flat 
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MAM 






A monopole-antimonopole chain with two monopoles 
on the symmetry axis, one above and one below the 
xy-plane, and an antimonopole at the origin. 


VMV 


c 




A monopole at the origin and two vortex rings, one 
above and one below the xy-plane. 


MV 




» : — o 


A monopole at the origin and a vortex ring in the 
xy-pl&ne. 


VMAMV 






Two vortex rings, one above and one below the 
xy-plane, two monopoles on the symmetry axis and 
an antimonopole at the origin. 


VMVV 






Three vortex rings, one above and one below the 
xy-plane and one in the xy-plane, and a monopole at 
the origin. 



Table 2: Configurations of the nodes and vortex rings of the Higgs field for m = 3, n = 3 
solutions. 



space solutions, for several values of A in this critical range. 

While the left part of the figure shows the evolution of the zeros of the Higgs field 
with increasing scalar coupling A, the right part gives for a particular value, A = 0.815, 
the assignment of the branches, to which the zeros belong. The a-branch emerging from 
the fundamental A-branch is labeled c, and the a-branches emerging from the new lower 
A-branch and the new upper A-branch are labeled a and b, respectively. 

It is intuitively clear that in the complicated pattern of interaction between the 
monopoles and antimonopoles, the gravitational attraction plays a role similar to the 
scalar interaction, thus one can expect that coupling to gravity may provide a com- 
pensation for the weakening of the scalar interaction. This is indeed, what we here 
observe regarding the existence of the solutions and their node structure. For increasing 
a, the branches a and b merge and disappear, while the branch c persists, but changes 
character, reversing the steps seen for the fundamental flat-space A-branch in Fig. [7J 

Indeed, at a = 0, all three branches begin as MAC solutions. As a increases, the two 
MAC branches a and b merge and end at a small a max - The branch c, however, extends 
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Figure 8: (a) The location of the isolated nodes z and vortex rings (p, z) of the Higgs field of 
m = 3, n = 3 solutions versus a for several values of the scalar coupling A in the critical range 
A Cl < A < A C2 and beyond, (b) The isolated nodes and vortex rings of the Higgs field for the 
complete set of solutions for A = 0.815. (See Table 2 for the notation of the node and vortex 
ring configurations of the Higgs field.) 

further and changes its node structure accordingly, as a increases. First a vortex ring 
emerges from each outer node. Next the nodes move towards each other, merging at the 
origin. Then a central vortex ring appears in the xy-plane, and finally the outer rings 
merge with the central ring. 

The branch c still persists to somewhat larger values of a, and then merges at an 
QWx with a second branch of vortex ring solutions labeled d. As expected, this branch 
extends backwards to a = 0, where it connects to a generalized BM solution. 

The evolution of the solutions with increasing A can now be understood from the left 
part of the figure. As A increases, the branches a and b extend to larger values of a, before 
they merge. At the same time the nodes of the branches b and c approach each other. 
At the critical value A C2 the two flat-space branches merge, and their nodes coincide. 
As A increases further, this critical point (where two branches merge) is retained, but it 
moves to finite values of a. Note, that the pattern of nodes of the gravitating solutions 
in the interval [0.03 < a < 0.07] for A = 0.8385 exhibited in Fig. [HJ nicely reflects the 
(reversed) pattern of nodes of the flat-space solutions in the interval [0.8 < A < 0.85] 
exhibited in Fig. [3 

Emergence of the new fiat-space branches at A C3 

As seen above, the bifurcations of the flat-space branches at A Cl and A C2 are reflected 
by bifurcations in curved space. In particular, the critical values of the gravitational 
coupling strength a, where bifurcations arise for a given value of the scalar coupling 
strength A, move towards smaller values of a with decreasing A. As the minimal critical 
value of a reaches zero, the flat-space bifurcation point A C2 arises. 

It is thus obvious, that the new flat-space A-branches arising at A C3 , should also have 
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(d) a (e) a 

Figure 9: The location of the isolated nodes z and vortex rings (p, z) of the Higgs field for the 
complete set of solutions for A = 1 (a) and A = 2 (b) versus a. Change of the structure of the 
nodes of the Higgs field following the evolution of the a-branches (c). The mass M (d) and 
the value of metric function at the origin 1(0) (e) for the same set of solutions. 
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precursors in curved space. These precursors indeed led to the detection of the new 
flat-space branches. In the following we analyze the branches of gravitating solutions, 
obtained with increasing scalar coupling strength A, to clarify the emergence of the new 
flat-space branches. 

Let us first consider the sets of solutions at A = 1 and A = 2, exhibited in Fig. [91 
While the node structure nicely reveals the types of solutions present at the respective 
value of A, and their evolution with a, the values of the metric functions at the origin are 
instructive to get an overview of the sets of solutions present. For better identification, 
we again label the various branches of solutions a, b, c, etc. 

For A = 1 there are four a-branches, just like for the smaller A values of Fig. [81 When 
moving continuously along these a-branches, the node structure changes from MAC, 
present for a and most of b, to MAC plus vortex rings, shortly before the bifurcation 
with c. The size and the location of the nodes and the rings then evolves along c until a 
single ring is left in the xy-plane (together with the always present node at the origin). 
This node structure is then retained on branch d, where the solutions evolve towards 
the generalized BM solution as a —>■ 0. 

As A increases further, interestingly, another bifurcation arises, and two more a- 
branches appear. Consequently, there is now an interval of a with six a-branches. The 
branches labeled a to f are exhibited in the figure for A = 2. Moving continuously along 
the branches, one again observes the structure of the solutions change from MACs to 
MACs plus vortex rings, and then to a monopole plus vortex ring(s). 

Most important for understanding the emergence of the new flat-space branches is 
the evolution of the branches d and e with A, with particular emphasis on the critical 
value of a, where these two branches merge. We therefore consider in Fig. [10] sets of 
solutions with still higher values of A. 



1=2.3 1=3.92 1=6.48 




0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 

(a) oc (b) a (Cj a 



Figure 10: The value of metric function at the origin, 1(0), for A = 2.3 (a) and A = 3.92 (b) 
and for A = 6.48 (c). 

As A increases further, the branches d and e extend to increasingly smaller values 
of a. At the critical value A C3 of the scalar coupling, the critical value of a, where the 
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branches d and e merge, then precisely reaches zero, giving rise to the bifurcation point 
of the new flat-space solutions. 

For values of the scalar coupling beyond A C3 , the set of solutions then splits into two 
disconnected parts. The first part emerges with branch a from the respective fundamen- 
tal flat-space solution, evolves via branch b and branch c, and reaches along branch d 
the new upper mass flat-space solution in the limit a — > 0. The second part emerges 
with branch e from the new upper mass flat-space solution, and evolves with branch f 
towards the limiting generalized BM solutions in the limit a — > 0. This is demonstrated 
in Fig. [10] for the set of solutions with A = 3.92. 

Solutions at large A 

In the case of the m = 2, n = 3 solutions we observed the appearance of new 
a-branches of solutions for large values of A, exhibiting an intriguing behaviour: In 
the limit A — ► oo the minimal value a m i n , beyond which these new branches exist, 
decreases towards zero, while the critical solutions themselves approach the generalized 
BM solution as a m j n — > oo. 

We here observe an analogous pattern for the m = 3, n = 3 solutions. To see this 
pattern arising, we reconsider the evolution of the branches a and b with A, exhibited 
in Fig. M and Fig. [TDJ As A increases, the branches a and b first elongate towards larger 
values of a. Then at a critical value of A, a new bifurcation arises, and two additional 
branches appear. The minimal value a m \ n possible for these new branches now decreases 
with increasing A, analogous to the case of the respective branches of the m = 2, n = 3 
solutions. 

3.4 Topologically trivial sector: m = 4, n = 3 
3.4.1 Flat space solutions 

We now turn to the m = 4, n = 3 solutions, which reside in the topologically trivial 
sector. As before we first address the flat-space solutions, reviewing and supplementing 
previous results [TO] . 

We exhibit in Fig. [JTj the mass and the node structure of the m = 4, n = 3 flat- 
space A-branches. The fundamental A-branch is present for arbitrary scalar coupling 
strength. Its solutions possess two vortex rings, located symmetrically with respect to 
the xy-plane. 

At A Cl = 0.491 a first bifurcation arises, where a pair of new A-branches appears, 
possessing higher mass and a mixed node structure, with two outer nodes on the sym- 
metry axis and a vortex ring in the xy-plane [TUJ. With increasing A the solutions on the 
lower mass branch retain their node structure, while the upper mass branch solutions 
transform first into MAC solutions with four isolated nodes on the symmetry axis, and 
then back to the mixed node structure. 



21 




MVA 



(a) ° 2 4 ^ 6 s 10 (b)0 e 

Figure 11: (a) The mass M of the flat-space solutions on the fundamental A-branch as well 
as on the additional A-branches for m = 4, n = 3 solutions, (b) The location of the isolated 
nodes z and vortex rings (p, z) of the Higgs field for the same set of solutions. (See Table 3 for 
the notation of the node and vortex ring configurations of the Higgs field.) 



MAMA 


: 




A chain with two monopoles and two 
antimonopoles in alternating order on the 
symmetry axis. 


MVA 






A monopole above and an antimonopole below 
the xy-plane on the symmetry axis and a 
vortex ring in the xy-plane. 


VV 






Two vortex rings, one above and one below the 
xy-plane. 


VV 






Two concentric vortex rings in the xy-plane. 



Table 3: Configurations of the nodes and vortex rings of the Higgs field for m = 4, n = 3 
solutions. 



Concerning the mass of these solutions, a transition between the fundamental branch 
and the new lower mass branch is observed at A^ 59.8, making the new lower branch 
solutions energetically favourable. As for the m = 2 and m = 3 solutions, it appears 



22 



energetically advantageous for the m = 4 solutions to exchange vortex rings for isolated 
nodes, when the scalar coupling is large. 

We therefore argued before |10] , that a further bifurcation was likely to exist, where 
another pair of branches would appear, representing MAC solutions with four isolated 
nodes on the symmetry axis. The solutions on this (conjectured) lower (mass) branch 
should then become the energetically most favourable configurations for large values of 
A. 

Indeed, as seen in Fig. [11], such a second bifurcation arises at A C2 = 1.288, and a pair 
of MAC solutions appears. With increasing A the solutions on the lower mass branch 
retain the MAC node structure, while the node structure of the upper mass branch 
solutions changes and forms two vortex rings. Concerning the mass of the solutions we 
also observe the conjectured transition: beyond A^ ~ 45.2 the new MAC branch has 
the lowest mass and thus becomes the energetically favoured branch. 

Surprisingly, at A C3 = 3.962 a third bifurcation arises and a third pair of branches 
appears with a new type of node structure. These solutions possess two concentric vortex 
rings in the xy-plane. Consequently their mass is considerably higher than the mass of 
the other configurations. 

But these two branches of solutions have another interesting feature. Unlike the 
other branches of solutions, this third pair of branches exists only in a relatively small 
range of the scalar coupling. At the bifurcation point A C4 = 66.354, the pair of branches 
merges and disappears again. 

Since this third pair of branches of solutions appeared unexpectedly, the existence 
of further branches of flat-space solutions in certain ranges of the scalar coupling seems 
possible. 

3.4.2 Gravitating solutions 

When gravity is coupled, branches of gravitating solutions arise from all of these m — 4, 
n = 3 flat-space configurations. But as seen above, such a-branches can bifurcate many 
times, leading to a plethora of gravitating solutions for larger values of A. We therefore 
refrain from obtaining the complete picture for these m = 4, n = 3 solutions, and 
present only results for a few selected values of A. In particular, we exhibit in Fig. [TJ] 
and Fig. [T5] the mass, the node structure and the metric function value 1(0) versus a, 
for sets of gravitating m = 4, n = 3 solutions with A = 1, 2 and A = 4.5. 

The first value A = 1 is chosen in the interval A Cl < A < A C2 , after the first bifurcation. 
So there are three flat-space solutions, each giving rise to an a-branch. The a-branch 
emerging from the fundamental flat-space solution connects via a second a-branch to 
the generalized BM solution, while the a-branches emerging from the first new pair of 
flat-space solutions merge with each other. By considering the evolution of the nodes 
on the new upper mass branch with A and with a, respectively, we note again, that an 
increase of the scalar coupling can lead to an analogous effect as the decrease of the 
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Figure 12: The mass M of m = 4, n = 3 solutions versus a for A = 1 (a), A = 2 (b). The 
location of the isolated nodes z and vortex rings (p, z) of the Higgs field for the same set of 
solutions for A = 1 (c), A = 2 (d), and likewise the value of the metric function at the origin 
1(0) for A = 1 (e), A = 2 (f). (See Table 3 for the notation of the node and vortex ring 
configurations of the Higgs field.) 
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gravitational coupling. In both cases, the MAC node structure can change into a mixed 
node structure at a critical value of the coupling. 

The second value A = 2 resides in the interval A C2 < A < A C3 , after the second 
bifurcation. Thus there are five flat-space solutions, each giving rise to an a-branch. As 
in the previous case, the a-branch emerging from the fundamental flat-space solution 
connects via a second a-branch to the generalized BM solution, and this feature of the 
fundamental branch solutions is retained, independent of A. 

However, the a-branches emerging from the first new pair no longer merge with each 
other. Instead the first new lower branch merges with the second new upper branch, 
and the first new upper branch merges with the second new lower branch. The reason 
for this can be inferred from the evolution of the a-branches with A, since the emergence 
of the second new pair of flat-space solutions can be traced back to the appearance of 
a bifurcation on one of the a-branches of the first new pair of solutions, for a value of 
A > 1. As A then increases, the two additional a-branches grow in size, while their 
bifurcation point a m i n moves towards smaller values of a, reaching a = at A C2 . Thus 
beyond A C2 , the a-branches of the first pair of solutions are disconnected from each other 
and merge instead with those of the second pair. 

As seen in the figure, at A = 2 the solutions on the a-branch emerging from the first 
lower mass solution possess a mixed node structure, those connected to the second upper 
mass solution have a node structure changing from MAC to mixed, while the solutions 
on the other branches have MAC structure. 

The third value of A is chosen in the interval A C3 < A < A C4 , where three pairs 
of flat-space branches are present. Their corresponding a-branches are exhibited in 
Fig. [131 At this value of A a plethora of a-branches is present, caused by numerous 
further bifurcations. The node structures of these solutions range from MACs to double 
vortex rings. Thus one may speculate about the existence of solutions with still more 
complicated node structures, which might arise from new a-branches, as the scalar 
coupling is further increased. 

4 Conclusions 

We have investigated static axially symmetric solutions of SU(2) Einstein- Yang-Mills- 
Higgs theory, representing monopole-antimonopole pairs, chains, vortex rings, and new 
types of configurations with mixed node structure. Such new configurations appear for 
larger values of the scalar coupling, when the subtle interplay between repulsive and 
attractive forces allows for more than one non-trivial equilibrium configuration of these 
systems. 

In flat space, at critical values of the scalar coupling, bifurcations arise, where pairs 
of new branches of solutions appear, which possess a different node structure than the 
solutions on the fundamental branches. In particular, new solutions appear, where 
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Figure 13: The mass M of m = 4, n = 3 solutions versus a for A = 4.5 (a). The value of 
metric function at the origin, 1(0), for the complete set of solutions related to the new flat space 
A-branches (b) and for the first (c), second (d) and third (e) pair of a-branches emerging from 
the flat space A-branches. The location of the isolated nodes z and vortex rings (p, z) of the 
Higgs field for the solutions emerging the first (d), second (e) and third (f) pair of a-branches. 
(See Table 3 for the notation of the node and vortex ring configurations of the Higgs field.) 
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vortex rings present in the fundamental solutions, are replaced by isolated nodes on the 
symmetry axis. For high values of A the new solutions with monopole-antimonopole 
chain structure have the lowest mass. 

While we have studied here in detail only the monopole-antimonopole systems with 
m = 2, 3, 4 and n = 3, we conjecture, that this phenomenon is not restricted to these 
particular systems but that it is of a more general nature, implicating an enormous 
richness of configuration space for high values of A and larger m and n. For m = 4 and 
n = 3, for instance, we observed already seven distinct flat-space solutions in a certain 
range of the scalar coupling. 

The coupling to gravity leads to additional attraction between the components of 
these monopole-antimonopole systems, shifting the balance of forces and thus affecting 
the possible equilibrium configurations. Whereas for small scalar coupling only two 
gravitating branches of solutions are associated with a flat-space solution, there is a 
surge of further branches when the scalar coupling becomes large. 

Particular of these branches can be shown to give rise to new pairs of flat-space so- 
lutions. Thus study of the Einstein- Yang-Mills-Higgs configuration space gives valuable 
insight into the configuration space of Yang-Mills-Higgs theory. Indeed, by exploiting 
this interrelationship, we were able to find the new flat-space solutions presented. 

But the Einstein- Yang-Mills-Higgs configurations are also related to Einstein- Yang- 
Mills solutions. First of all, in the limit a — > 0, one of the a-branches of a particular 
set of solutions, characterized by m, n and A, always connects to a generalized Bartnik- 
McKinnon solution [TjJJ [201 EE] • But another connection arises in the limit of infinite 
scalar coupling. In this limit the bifurcation point a m ; n of two of the a-branches is 
seen to approach a = 0, while the solutions at a m i n approach the respective generalized 
Bartnik-McKinnon solution. 
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